Extracting from small lattices: unquenched results 
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We calculate the response of the microscopic Dirac spectrum to an imaginary isospin chemical 
potential for QCD with two dynamical flavors in the chiral limit. This extends our previous calcu- 
lation from the quenched to the unquenched theory. The resulting spectral correlation function in 
the e-regime provides here, too, a new and efficient way to measure Ftt on the lattice. We test the 
method in a hybrid Monte Carlo simulation of the theory with two staggered quarks. 
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I. INTRODUCTION 

A primary challenge for lattice gauge theory is the 
numerical determination of the low-energy constants of 
QCD with nearly masslcss u and d quarks. The chi- 
ral limit is a notoriously difficult problem on the lat- 
tice. Furthermore, one also has to face the issue of finite- 
volume effects, a problem that is compounded by the 
light quarks. A situation like this is best tackled head-on 
by means of finite-size scaling techniques, as is routinely 
done when the correlation length becomes comparable to 
the system size (such as near continuous phase transi- 
tions). Unfortunately, exact finite-size scaling functions 
for the chiral limit are not easy to come by. One notable 
exception is the so-called e-regime of QCD f ] . ft has been 
known for more than a decade that this makes pos- 
sible precise measurements of the infinite-volume chiral 
condensate S — {tpi-') from simulations on small lattices. 
By definition, in the e-regime the linear extent L of the 
box is smaller than the Compton wavelength 1 /m.^^ of the 
pions, but much larger than 1/Aqcd- It is in this sense 
that the lattice required is small: Rather than demand- 
ing that the pion be contained by the lattice, one requires 
precisely that the lattice be too small for the pion. 

In the e-regime the leading term of the effective par- 
tition function is a known function that is determined 
uniquely by the flavor symmetries and their spontaneous 
breaking. While this leads to an equally computable an- 
alytical expression for the condensate {'4>'4') that can be 
compared to lattice measurements 0] , a far more precise 
estimate comes from extracting S from the analytical 
expressions for the distributions and correlations of the 
smallest eigenvalues of the Dirac operator. The univer- 
sal predictions for different eigenvalue distributions ac- 
cording to the topological index v and gauge group rep- 
resentation of the quarks |3| have been checked in de- 
tail for both staggered and overlap fermions in a variety 
of studies with different color representations, mainly in 
quenched calculations (see, e. g., Refs. yJll for a partial 
list), but also with dynamical fermions [a,l3- 

Here we will consider the computation of the low- 



energy constant Ft^, associated with the decay of pi- 
ons. In the conventional approach _FV is measured in 
a regime where the Compton wavelength of the pion is 
much smaller than the lattice size (see 8] for some recent 
quenched and j^J for some recent dynamical calculations 
of F.^). We have recently presented a new method for ex- 
tracting the pion decay constant from Dirac operator 
eigenvalues in the e-regime [lOj |. In our first paper the 
focus was on quenched QCD. In this article, we extend 
this method to QCD with two dynamical light flavors, 
the physical u and d quarks. 

In principle, all low-energy constants of QCD can be 
determined numerically from the spectrum of the Dirac 
operator. To leading order the spectrum depends only 
on the infinite-volume chiral condensate E, but higher 
order terms (in the counting of the e-expansion for the 
chiral Lagrangian U) give rise to measurable deviations. 
The idea of the improved method proposed in Ref. [loj 
is to make use of an enlarged set of couplings provided 
by external fields. In particular we introduce a constant 
Abelian field that couples with opposite sign to the u and 
d quarks. This can be interpreted as an imaginary isospin 
chemical potential [Ol . Such a source turns out to have a 
dramatic effect on a certain spectral correlation function 
of the Dirac operator. It is this effect which allows for a 
very precise determination of F-^. 

An alternative method is to use a real baryon chemi- 
cal potential to extract the i^7r-dependence of the micro- 
scopic spectral density. Recently quenched results based 
on this method have been reported The a naly tical 
predictions for this method were derived in Ref. [ij, 13 ■ 
The advantage of our method is that the Dirac operator 
eigenvalues remain real, while a baryon chemical poten- 
tial makes them (and the determinant) complex. The 
pion decay constant has also been measured from meson 
correlation functions in the e-regime (see e.g. [iS]). 

The outline of the paper is as follows. We first present 
the analytical result for the spectral correlation function 
whose dependence on Ft^ is the backbone of our method. 
We then compare this prediction to the results of numer- 
ical simulations of QCD with two dynamical flavors, and 
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use this to extract the value for F-^. 



II. TWO LIGHT FLAVORS: THE ANALYTICAL 
RESULT 



otherwise artificial use of an imaginary isospin chemical 
potential. Viewed in terms of external gauge potentials 
it is actually the natural choice, corresponding to a con- 
stant real Abelian gauge potential = Miso- 



An imaginary isospin chemical potential leads us to 
consider two distinct Dirac operators and their eigenval- 
ues for a given gauge potential A: 

^ (0(A)+zM)A"^-*A^"Vi"\ (1) 

i?_VA"' ^ (0(A) ~ z/.i,o7o)^i"^ = *AL"Vi"^ (2) 

Since the operators D± are anti-hermitian, the eigenval- 
ues A^''' and a!^"* lie on the real axis. This justifies the 



We consider the mixed two-point spectral correlation 
function of the Dirac operators D± . The formulae for the 
quenched theory were given in Ref. To]; here we briefly 
present the corresponding expressions for two dynamical 
flavors. The mixed two-point function, which depends on 
the two light masses to„ and and on /iiso, is defined 
as 



n I 71 I 

where the brackets denote the average over the QCD partition function with two flavors and fiiso 7^ 0, 
^2(to„, mrf; i/iiso) = / [d^]^ det(D+ + m„) det(i:'_ -I- ma) e-^™^'^). 



(3) 



(4) 



The index on the measure indicates that this is the partition function in a sector where the gauge fields have 
topological charge i^. Keeping fixed the scaling variables 



as 1^ — *■ 00, we define the microscopic two point correlation function 



(5) 



pI^^ i^-,'rhu,rnd] i/tiso) = Hni 



_! (2)( ^ ^ ■ Ai: 



E2y2 ^ \j:v'y.v'y.v'y.v' F^W 



(6) 



This spectral correlation function p'^) q^j^ obtained from the mixed "partially quenched" scalar susceptibility. 



xim+,m^,mu,md;ifiiBo) = ( tr 
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which possesses the spectral representation 



n «A^ ' + 



I iX'^S' + TO- 



The desired correlation function follows from the double discontinuity across the imaginary axis |l6j | , 

1 



p(2)(A+,A_ 

; ^d-i ^Miso j 



rDisc x(to-i-, to_, to„, to^; i^iso)| ™+='A+ 



(7) 



(8) 



(9) 



47r2 
1 

-x(«A+ + e, - e, to„, irid] ifJ-iso) + x(jA+ - e, - e, to„, to^; iMiso)]- (10) 
I 



2 lim [x(iA+ + e, iA_ + e, to„, to^; i^iiso) - x(«A+ - e, iA_ + e, to„, to^; i/Ziso) 



The power of the e-regime is that this function, p^^\ can be expressed in closed analytical form. 
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The mixed scalar susceptibility lO involves two addi- 
tional fermion species, with mass m_|_ and m_, coupled 
as well to the external field /iiso- This quantity thus has 
to be evaluated in a theory with four fermion species, of 
which only two become physical (the u and d quarks), 



while the remaining two are quenched away and serve as 
sources for the two eigenvalues involved in the two-point 
correlation function. In order to do this we will use the 
replica method (see, e.g., 0|). We have 



x(m+,TO_,m„,md;i//iso) = lim ^9,„^5„i_ log Z2„^2("i+, to_, m„, m^; i//iso), 



n-+0 n 



(11) 



where the generating functionals 2^2n,2 have 2n replica flavors besides the two flavors of mass m„ and m^, 

X det(I?+ + m„) det(i:>_ + m^) e"^^"*^^ 



(12) 



Half of the replica flavors have mass m+ and chemical potential +i/iiso while the other half have mass m_ and chemical 
potential —i^iso- 

In the e-regime the leading terms of the effective partition functions satisfy a string of Toda lattice equations which 
can be used to extract the replica limit 0, 0| . Specifically [31 j 



rh+drn+m_drh_ log Z2„,2(to+, TO-, ?7ld; ifliso) 



4n^(?7t+?7i_) 



2 ^2n+2,2 



{m+, m_, m„, m^; «Aiso)^2n-2,2(™+, W-, rhd, «Aiso) 



[Z2„,2("i+, m-, m„, rfid; iftisoW 
Taking the replica limit (|ll|l of this Toda lattice equation we get 

•^2,2(™+, TO-, m«, rhd.; ifiiso)Z2-2{mu, TOd|m+, 7fi_; i/iiso) 



Xim+,m-,mu,rhd-, ifliso) = 4m+m_ 



[Z2{mu,md]ifiisoW 



(13) 



(14) 



We need three different effective partition functions to evaluate the right hand side of this equation. Of these, the 
two that involve only positive numbers of flavors are of the form 



^2n,2(TO+,TO-,TO„,md;zAiso) 

IUeUi2n+2) 

where B = diag(l„+i, — l„+i) and AI = diag(m 



dL/det(C/)''e3^^'^-°'^''[^'-^l['^*'-^]+^^'^'^'(*^^'^+*^'^^^ 
m^,md). In particular, 



Z2{m+,m^; ifliso) = e 



dtte^^"°^ h{trh+)h{trh- 



(15) 



(16) 







where we have defined rh± — m±TiV , while an explicit analytical expression for ^2,2 can be obtained from Ref. |20j, 
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^2,2(to+, TO_, m„, mrf; ifi-^so) 



(m^ — m^)(m_ — m^) 



Z2(m+,m_;i/iiso) 2'2(to+, m^; «Aiso) 
Z2{mu,'m^; ifliso) ^2( 



(17) 



The remaining partition function Z2-2 on the right hand side of Eq. p4|) is the most difficult part. It involves two 
fermion species (the two physical u and d quarks) and in addition —2 flavors, i.e., two bosons. The resulting effective 
axtition function is an integral over a graded group. Deferring the technical discussion of this to a forthcoming paper 
we quote the final result. 



•^2,-2(m„, md|m+, i/iiso) 



(m^ - m\)Z2{mu, rhd; ifliso) Zi _i(md|TO_) 

-2'i._i(m„|m+) {ml - m^)Z_2(TO+, W-; «/tiso) 



In Eq. H18|l the (/^-independent) partition function with one fermion and one boson is given by 

Zi-i{x\y) = xIi{x)KQ{y) + yKi{y)Io{x), 



(18) 



(19) 
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and the partition function with two bosons is 



(20) 



where x = xT,V and y = yTiV. 

We can finally compute the required two-point correlation function © by means of the discontinuity (PJ. The 
result is 



(21) 



[m„/^+i(?7i„)J,,(^+) + ^+J^+i(^+)/y(7fi„)] ['mdI„+i{md)J^{(.-) + £.-J^+ii£.-)Limd)] 



For a numerical evaluation it is advantageous to rewrite the improper integral appearing in the fourth line as 

/ dUe-2ALtV,(te+)J,(<e-) = T^exp (-^±±^\ (t^) - I dtte-^^^"^^'" Mt^+)Mt^^). (22) 



In exact analogy to the quenched case 0|, the corre- 
lation function H21() changes dramatically when /ijso is 
made non-zero. A fit to Monte Carlo data with /ijso 7^ 
using Eq. H21|l will then readily produce a measurement 
of iv- 



III. NUMERICAL RESULTS 

We demonstrate our method by applying it to the 
eigenvalue correlation function calculated for SU(3) lat- 
tice gauge theory with two staggered fermions of equal 
mass m„ = uid = m. The gauge action is the standard 
plaquette action and unimproved staggered operators are 
used throughout. This theory has an exact U(2)xU(2) 
symmetry; while it represents eight tastes in the contin- 
uum limit, we work far from the continuum limit, at a 
strong coupHng f3 = 4.2, in order to stay on the confin- 
ing side of bulk and finite-temperature phase transitions 
|22j |. In addition, and since we only aim at illustrating the 
method here, this rather strong coupling provides us with 
an unambiguous number of pseudo-Goldstone modes. At 
the coupling considered, two staggered fermions give rise 
to spontaneous chiral symmetry breaking of the form 
U(2)xU(2)— »U(2). The staggered theory at this couphng 
thus appears to be "insensitive" to gauge field topology 
(see for instance jl^): The additional U(l) factor in the 
coset dictates that we must compare our numerical re- 
sults with = in Eq. (|^ . 



We ran the exact Hybrid Monte Carlo algorithm |2^ on 
lattices with volume F = 4^ and 6^, setting /iiso = 0.0125 
and 0.0055, respectively, in order to keep constant the 
scaling variable /tjso- We studied the cases m — 0.01 
and 0.025 for the smaller volume, and correspondingly 
m — 0.002 and 0.005 in the larger volume. Note that the 
scaling variable rh = mT,V is fixed and mT^L ^ 1 (cf. ta- 
ble^ . An idea of the size of the simulation can be gotten 
from the run parameters listed in Table We note that 
our method requires very small quark masses and hence 
many CG iterations (typically 10"* for the smaller mass in 
V — 6'^, restarting every 500) for inverting the fermion 
matrix; the restriction to small volume, while allowing 
easy replication of the simulation for many processors, 
precludes gaining any advantage from true parallel com- 
putation. 

As discussed in when extracting physical observ- 
ables we must keep in mind that the continuum theory 
describes 4 tastes of quark. Therefore, to determine the 
values of and S from the staggered eigenvalue spec- 
trum we must replace V in the analytical predictions by 
W. 

We begin with our main numerical result, the pion 
decay constant F^^. In order to improve the statistics we 
consider the integrated correlation function 

Pint (2;) EE y d^p'^^\^ + x,^,mu,md-,ifi.iso)- (23) 
In Fig. n we show this integrated correlation function 
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V m trajectories dt steps/trajectory trajectories/measurement 

■7 aOl 4000 0067 l50 i 

4^ 0.025 4000 0.0125 80 1 

6'^ 0.002 27960 0.001 100 5 

e** 0.005 5135 0.001 200 5 



TABLE I; Run parameters for the Hybrid Monte Carlo simulations. 




FIG. 1: The integrated correlation function 12311 with ^min ~ 5 and ^max = 55 for the ensembles with 1/ = 6* and /iiso = 0.0055. 
Left: m = 0.002. The curve is the resuh of a fit with E = 0.3200 ± 0.0005 and Fn = 0.3379 ± 0.0018, with x'/dof = 1.33. 
Right: m = 0.005. The fit parameters are E = 0.3212 ± 0.0012 and F.^ = 0.3384 ± 0.0025, with ^ j do] = 1.13. 



as measured on our 6"* ensemble, for both masses. The 
curves are fits to Eq. via Eq. ||2U , see Table |nl 
The weighted averages of the fit parameters for the two 
data sets are S = 0.3204 ± 0.0004 and = 0.3381 ± 
0.0010, in lattice units. There are a number of other 
ways to determine S from the eigenvalue densities and 
correlations; we will see that the value just obtained is 
consistent with other features of the spectra. 

We now proceed to study the sensitivity of the correla- 
tion function to the mass. Since the mass affects mainly 
low-lying eigenvalues, this sensitivity is greatest when the 
correlation function is not integrated. We show in Fig. [21 
the measured correlation function with one eigenvalue 
fixed at ^_ = 4. For each data set we compare to the- 
oretical curves [Eq. if^ ] with masses corresponding to 
the two chosen values of m in the simulation, as well as 
to infinitely high mass (the quenched theory). We set S 
and Ft^ to the values obtained in the fits shown in Fig. ^ 
for each mass separately. It is clear that the effect of 
the mass in the analytical prediction is consistent with 
the data and inconsistent with the quenched result. It 
would take, however, a somewhat greater number of con- 
figurations to get a precise test of the mass dependence. 

The scaling properties of the e regime can be studied 



by comparing distributions for the two volumes. We be- 
gin with the eigenvalue density shown in Fig. |3| For the 
larger volume (right hand figure) the density as expected 
rises to a constant level (by convention the value reached 
should be I/tt as indicated by the horizontal line). The 
fall-off around the value ^ = 80 occurs since only the 
lowest 24 eigenvalues have been computed. On the other 
hand the eigenvalue density in the smaller volume (left 
hand figure) only briefly remains at the plateau of height 
I/tt. Thus it is not advisable to consider the integrated 
correlation function over the entire range of the 24 eigen- 
values computed. In Fig. 0] we consider rather the cor- 
relation function with ^_ fixed at the value 4. Since we 
have chosen the quark masses such that mV is fixed, the 
analytical curves shown are identical to those in Fig. 12 
The general agreement confirms the expected scaling be- 
havior. 

Finally we take a brief look at the case where /i — 
0. In this case the operators Dj^ and D_ are identical 
and the correlation function consequently has a delta- 
function peak at equal arguments. In Fig. [3 we show the 
integrated correlation function at /ijso = as obtained on 
the 6* lattice at mass m = 0.002 along with the analytical 
prediction, using the value of S determined in Fig.^ The 
consistency with the data is a non-trivial cross check of 
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V 


m 




E 


m,r / Ftt 






0.002 


0.3379 ±0.0018 


0.3200 ± 0.0005 


0.3133 


0.6354 


6* 


0.005 


0.3384 ± 0.0025 


0.3212 ±0.0012 


0.4949 


1.005 



TABLE II: Results of the data analysis. To extract the ratio m^/F^ we used the GOR-relation m^F^ — 2mE. 




0246802468 

FIG. 2: The correlation function Q, with ^„ fixed at 4, measured in the ensembles with V — 6*. Left: m = 0.002. Right: 
m = 0.005. The curves are Eq. ^ with m = 0.002 (solid), 0.005 (dashed), and oo (dotted). 




10 20 30 40 50 60 70 80 90 10 20 30 40 50 60 70 80 90 

FIG. 3: The eigenvalue density of the 24 lowest eigenvalues measured for the smaller scaled mass: V = 4*, m = 0.01 (left) and 
V = 6'*, m = 0.002 (right). The horizontal line in each case is at ps — I/tt. For 1/ = 4* the eigenvalue density only remains 
briefly at the height l/vr whereas for V — 6* all of the first 24 eigenvalues are within this plateau. 



the value of E obtained in Fig. The delta-function at 
^+ = ^_ is not shown. As /iiso is turned to a non-zero 
value the deha function peak gets smeared into a broader 
peak around zero as in Fig. ^ It is the fit to the shape 
of this broader peak that allows us to determine F^. 



IV. SUMMARY 



We have demonstrated that F^^ can be extracted with 
high precision from small lattices in unquenched lattice 
QCD. The approach exploits the dramatic dependence 
of the universal microscopic correlation functions on F.^. 
This dependence enters through the coupling to an ex- 
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FIG. 4: The correlation function I^J, with ^_ fixed at 4, for V = i'^; m = O.Of (left) to be compared with the full curve and 
m = 0.025 (right) to be compared with the dashed curve. The curves are obtained from 1211 with E and -FV as determined in 
Fig.0 The dotted curve indicates the quenched result. 




has allowed us to extract F^r in lattice units. See Table 
^for a summary. 

In this first study of the method with dynamical 
fermions we have worked with an unimproved staggered 
action. It would be of great interest to perform a sim- 
ulation with actions suffering less severely from scaling 
violations, thus allowing to quote a result for in phys- 
ical units. 
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FIG. 5: The integrated correlation function at fUso = 0, for 
V — m = 0.002. The delta function at i = is not shown. 
The comparison with the analytic curve uses the value of E 
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ternal vector field in form of an imaginary isospin chem- 
ical potential. The derivation of the analytical results 
used here is quite involved and will be presented in [2H . 
We have performed appropriate two flavor lattice sim- 
ulations, and the numerical results have been shown to 
agree with our analytical predictions. This agreement 



PHD would like to thank the Yukawa Institute for The- 
oretical Physics at Kyoto University where this work was 
completed during the YITP workshop on "Actions and 
Symmetries in Lattice Gauge Theory" YITP-W-05-25. 
The work of KS we financed by the Carlsberg Foun- 
dation. The work of BS was supported in part by the 
Israel Science Foundation under grant no. 173/05. He 
thanks the Niels Bohr Institute for its hospitality. The 
work of DT was supported by the NSF under grant 
No. nsf-phy0304252 He thanks the Particle Physics 
Group at the Rensselaer Polytechnic Institute for its hos- 
pitality. Our computer code is based on the public lattice 
gauge theory code of the MILC Collaboration [25| . 



[1] J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987); 
188, 477 (1987); H. Neuberger, Phys. Rev. Lett. 60, 889 
(1988); H. Leutwyler and A. Smilga, Phys. Rev. D 46, 



5607 (1992). 

[2] E. V. Shuryak and J. J. M. Verbaar schot, Nucl. Phys. 
A 560, 306 (1993) |hep-th/9212088l ; M. E. Berbenni- 



8 



Bitsch et ai, Nucl. Phys. Proc. Suppl. 63, 820 (1998) 
[hep-lat/9709102 ; P. H. Damgaard, U. M. Heller 
and A. Krasnitz, Phys. Lett. B 445, 366 (1999) 
lhep-la t/9810060 . 

[3] J. J. M. Verbaarschot, Phys. Lett. B 368, 137 (1996) 
piep-ph/9509369 ; P. H. Damgaard, R. G. Edwards, 
U. M. Heller and R. Narayanan, Phys. Rev. D 61, 094503 
(2000) hep-lat/9907016 ; P. Hernandez, K. Jansen 
and L . Lellouch, Phys. Lett. B 469, 198 (1999) 
liiepT at/9907022 . 

[4] P. H. Damgaard and S. M. Nishigaki, Phys. Rev. D 63, 
045012 (2001) hep-th/0006111 . 

[5] M. E. Berbenni-Bitsch, S. Meyer, A. Schafer, 
J. J. M. Verbaarsc hot and T. We ttig, Phys. Rev. 
Lett. 80 (1998) 1146 |hep-lat/97040T8| ; R. G. Edwards, 
U. M. Heller, J. E. Kiskis and R. Narayanan, ibid. 
82, 4188 (1999) hep-th/9902117l; R. G. Edwards, 
U. M. Heller and R. Narayanan, Phys. Rev. D 60, 
077502 (1999) hep-lat/9902021 ; P. H. Damgaard, 
U. M. Heller, R. Niclasen and B. Svetitsky, Nucl. 
Phys. B 633, 97 (2002) hep-lat/0110028 |; L. Giusti, 
M. Liischer, P. Weisz and H. Wittig, J . High Energy 
Phys. 0311, 023 (2003) hep-lat/0309189'. 

[6] M. E. Berbenni-Bitsch, S. Meyer and T. Wettig, Phys. 
Rev. D 58, 071502 (1998) hep-lat/9804030 . 

[7] P. H. Damgaard, U. M. Heller, R. Niclasen and 
K. Rummukainen, Phys. Lett. B 495, 263 (2000) 
hep-lat/0007041 . 

[8] S. Aoki et ai, Phys. Rev. D 62, 094501 (2000) 
hep-lat/9912007 ; S. J. Dong, et al., Phys. Rev. 
D 65, 054507 (2002) hep-lat/0108020 ; T. W. Chin 
and T. H. Hsieh, Nucl. Phys. B 673, 217 (2003) 
hep-lat/0305016 ; S. Diirr and C. Hoelbling, Phys. Rev. 
D 72, 071501 (2005) hep-ph/0508085 ; C. Gattringer, 
P. Huber and C. B. Lang,Phys. Rev. D 72, 094510 (2005) 
hep-lat/0509003 . 

[9] C. R. AUton et ai, Phys. Rev. D 70, 014501 
(2004) hep-lat/0403007 ; F. Farchioni, L Montvay 
and E. Scholz, Eur. Phys. J. C 37, 197 (2004) 
hep-lat/0403014 ; C. Aubin et ai, Phys. Rev. D 
70, 114501 (2004) hep-lat/0407028 ; M. F. Lin, PoS 
LAT2005, 094 (2005) hep-lat/0509178 . 
[10] P. H. Damgaard, U. M. Heller, K. SplittorfT and B. Svetit- 
sky, Phys. Rev. D 72, 091501 (2005) hep-lat/0508029 . 
[11] T. Mehen and B. C. Tiburzi, Phys. Rev. D 72, 014501 



(2005) hep-lat/0505014'. 
[12] J. C. Osborn and T. Wettig, Proc. Sci. LAT2005, 200 

(2005) hep-lat/0510115 . 
[13] G. Akemann and G. Vernizzi, Nucl. Phys. B 660 (2003) 

532 hep-th/0212051 ; J. C. Osborn, Phys. Rev. Lett. 93, 

222001 (2004) hep-th/0403131 ; G. Akemann, J. C. Os- 
born, K. SplittorfT and J. J. M. Verbaarschot, Nucl. Phys. 

B 712 (2005) 287 hep-th/0411030 . 
[14] K. SplittorfT and J.J.M. Verbaarschot, Nucl. Phys. B 683, 

467 (2004), hep-th/0310271 . 
[15] W. Bietenholz et al., JHEP 0402, 023 (2004) 
hep-lat/03 11012 ; L. Giusti, P. Hernandez, M. Laine, 

P. Weisz and H. Wittig, JHEP 0404, 013 (2004); 

jliep-lat/0402002 . H. Fukaya, S. Hashimoto and 

K. Ogawa, Prog. Theor. Phys. 114, 451 (2005) 

"hep-lat/0504018 ; S. Shcheredin and W. Bietenholz, PoS 

LAT2005, 134 (2005) hep-lat/0508034 . 
[16] P. H. Damgaard, J. C. Osborn, D. Toublan and 

J. J. M. Verbaarschot, Nucl. Phys. B 547, 305 (1999) 
hep-th/9811212 ; D. Toublan and J. J. M. Verbaarschot, 

tbid. 603 (2001) 343 hep-th/0012144 . 
[17] P. H. Damgaard and K. SplittorfT, Phys. Rev. D 62, 

054509 (2000) hep-lat/0003017 . 
[18] E. Kanzieper, Phys. Rev. Lett. 89, 250201 (2002) 

cond-mat/0207745 . 
[19] K. SplittorfT and J.J.M. Verbaarschot, Phys. Rev. Lett. 

90, 041601 (2003) cond-mat/0209594 ; Nucl. Phys. B 

695, 84 (2004) hep-th/0402177 . 
[20] G. Akemann, Y. V. Fyodorov and G. Vernizzi, Nucl. 

Phys. B 694 (2004) 59 hep-th/0404063 . 
[21] P. H. Damgaard, U. M. Heller, K. SplittorfT, B. Svetitsky, 

and D. Toublan, in preparation. 
[22] J. B. Kogut and D. K. Sinclair, Nucl. Phys. B 295, 465 

(1988); F. R. Brown, et al, Phys. Rev. D 46, 5655 (1992) 
hep-lat/9206001 . 
[23] P. H. Damgaard, U. M. Heller, R. Niclasen and 

K. Rummukainen, Phys. Rev. D 61, 014501 (2000) 
hep-lat/9907019 . 
[24] S. Duane, A. D. Kennedy, B. J. Pendleton and 

D. Roweth, Phys. Lett. B 195, 216 (1987); R. Gupta, 

G. W. Kilcup and S. R. Sharpe, Phys. Rev. D 38, 1278 

(1988). 
[25] Available Trom 

[http : / /www ■ physics . Utah . edu/~ detar/milc/| 



